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RECURSIVE KALMAN FILTER FOR 
FEEDBACK FLOW CONTROL IN REM FUEL CELL 

FIELD OF THE INVENTION 
[0001] The present invention relates to fuel cells, and more particularly 
to flow control in a fuel cell system. 

BACKGROUND OF THE INVENTION 

[0002] Fuel cell systems include a fuel cell stack that produces 
electrical energy based on a reaction between a hydrogen-based feed gas (e.g., 
pure hydrogen or a hydrogen reformate) and an oxidant feed gas (e.g., pure 
oxygen or oxygen-containing air). The hydrogen-based feed gas and oxidant 
feed gas are supplied to the fuel cell stack at appropriate operating conditions 
(i.e., temperature and pressure) for reacting therein. The proper conditioning of 
the feed gases is achieved by other components of the fuel cell stack to provide 
the proper operating conditions. 

[0003] The fuel cell system includes a compressor for compressing the 
oxidant feed gas to an appropriate operating pressure for reaction in the fuel cell 
stack. The compressor also supplies the oxidant feed gas to the fuel cell stack at 
a flow rate, which is dependent on the electrical load demand from the fuel cell 
stack. For example, when electrical load demand is increased, the compressor 
supplies the oxidant feed gas to the fuel cell stack at a higher flow rate. 

[0004] A mass flow meter is used as a feedback sensor for compressor 
mass flow control. The flow sensor monitors the mass flow rate of oxidant feed 
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gas through a pipe to the fuel cell stack. Typical mass flow meters are of a hot- 
wire type that are effected by the flow therethrough. Laminar flow through the 
pipe results in a steady, accurate signal from the mass flow meter. Turbulent 
flow, however, results in highly-varying, inaccurate signals. As a result, up to ten 
pipe diameters of straight pipe length is required before and after the mass flow 
meter to ensure laminar flow through the mass flow meter. In most applications, 
however, there is insufficient space to include such long sections of pipe. 

[0005] One conventional solution includes processing the flow sensor 
signal through a low pass filter. In this manner, the mass flow meter signal is 
smoothed, eliminating some of the erroneous data. This, however, still results in 
an inaccurate mass flow meter signal. Another conventional solution includes 
incorporating flow straighteners, such as honeycomb flow straighteners, within 
the pipe. As a result the component count is increased making the overall 
system more complex and expensive. 

SUMMARY OF THE INVENTION 
[0006] Accordingly, the present invention provides a method of 
operating a compressor of a fuel cell system. The method includes modeling a 
flow meter that measures a mass flow from a compressor with a first 
mathematical formula, generating a measured signal from the flow meter. The 
first mathematical formula and the measured signal are processed through a KF- 
based signal processing algorithm to provide a future signal estimate. The 
compressor is operated based on the future signal estimate. 
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[0007] In one feature, the first mathematical formula consists of a 3^ 
order model of the cathode system. 

[0008] In another feature, the method further includes predicting a 
current signal estimate based on a previously smoothed signal estimate. The 
previously smoothed signal estimate is determined based on a previously 
predicted estimate, a previous signal measurement and a previous gain. 

[0009] In yet another feature, the method further provides calculating a 
smoothed current signal estimate based on a predicted current estimate, a 
current measurement and a gain, the future signal estimate is based on the 
smoothed current signal estimate. The future signal estimate is further based on 
a current command signal. The method further includes modeling a compressor 
command signal with a second mathematical formula and calculating the current 
command signal based on the second mathematical formula. 

[0010] The present invention also include a fuel cell signal which is 
adapted to employ the above described methods for operating a compressor of a 
fuel cell. 

[0011] Further areas of applicability of the present invention will 
become apparent from the detailed description provided hereinafter. It should be 
understood that the detailed description and specific examples, while indicating 
the preferred embodiment of the invention, are intended for purposes of 
illustration only and are not intended to limit the scope of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0012] The present invention will become more fully understood from 

the detailed description and the accompanying drawings, wherein: 

[0013] Figure 1 is a schematic illustration of a cathode loop of a fuel 

cell system; 

[0014] Figures 2A and 2B are a set of graphs illustrating a compressor 
command signal generated by a controller of the fuel cell system and an air mass 
flow signal generated by a mass flow sensor of the fuel cell system; 

[0015] Figures 3A-3C are a set of graphs illustrating the compressor 
command signal and mass flow data points including a model overlay; 

[0016] Figure 4 is a flowchart illustrating the Kalman Filter (KF) -based 
mass flow signal processing according to the present invention; 

[0017] Figures 5 A and 58 are a set of graphs illustrating the 
performance of a recursive Kalman filter in accordance with the present 
invention; 

[0018] Figures 6A and 68 are a set of graphs comparing the KF-based 
mass flow signal processing and a traditional low pass mass flow signal 
processing; 

[0019] Figures 7A-7C are a set of graphs illustrating a compressor 
command signal and a mass flow signal generated using the KF-based signal 
processing of the present invention; and 

[0020] Figure 8 is a software diagram for the implementation of a 
recursive Kalman filter in accordance with the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0021] The following description of the preferred embodiments is 
merely exemplary in nature and is in no way intended to limit the invention, its 
application, or uses. 

[0022] Referring now to Figure 1 , a cathode loop of a fuel cell system 
10 is schemaitically illustrated. The fuel cell system 10 includes a fuel cell stack 
12 that processes hydrogen and oxygen to produce electrical energy. More 
specifically, as hydrogen-containing feed gas flows into an anode side of the fuel 
cell stack 12. A catalyst facilitates separation of the hydrogen-containing feed 
gas into electrons and hydrogen ions (i.e., protons). The hydrogen ions pass 
through an electrolyte membrane and combine with oxygen in a cathode side to 
produce water (H2O). The electrons, which cannot pass through the electrolyte 
membrane, flow from the anode side to the cathode side through an external 
circuit (e.g., the load). The load consumes the power generated by the fuel cell 
stack 12. 

[0023] A control element such as compressor 14 provides a feed 
stream such as oxygen-containing air for processing in the fuel cell stack 12. A 
mass flow sensor 16 measures the mass flow of air to the fuel cell stack 12. A 
controller 18 controls operation of the fuel cell system 10 based on control 
algorithms and an input 20. The input 20 generates a control signal based on a 
desired power output from the fuel cell system 10. For example, in the case of a 
vehicle, the input 20 is a throttle position signal and the control signal indicates a 
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desired vehicle speed. The controller 18 operates the fuel cell system 10 to 
obtain the desired power output from the fuel cell stack 12 for achieving the 
desired vehicle speed. The controller 18 controls the compressor 14 based on 
the feedback signal generated by the mass flow sensor 1 6. As described herein, 
the present invention is used to operate and control a compressor in the catalyst 
feed stream. However, one skilled in the art will recognize that the present 
invention could have utility for operating or regulation of other control alerts or 
measurements within the fuel cell system 

[0024] A set of graphs respectively illustrate an exemplary compressor 
command signal generated by the controller 18 in Figure 2A and an exemplary 
mass flow signal generated by the mass flow sensor 16 in Figure 2B. Initially, 
the compressor command signal is approximately 1,000 RPM. This corresponds 
to approximately 2 grams per second (gps) mass flow of air to the fuel cell stack 
12. At such a low mass flow rate, the air flow is laminar generating only slight 
noise in the sensor signal. After approximately 90 seconds of operating time, the 
controller 18 increases the compressor command signal to approximately 10,000 
RPM in step-wise fashion, corresponding to approximately 23 gps mass flow of 
air to the fuel cell stack 12. Such an increase can occur based on a desired 
power output increase from the fuel cell stack 12. The increased mass flow rate 
results in turbulent air flow which generates significant noise in the sensor signal. 

[0025] For example, at the lower air mass flow rate of about 2 gps the 
noise corruption is much smaller, with a variance of about 0.01 15, as compared 
to the higher air mass flow rate of about 23 gps, with a variance of about 3.636. 



6 



Attorney Docket No. 8540G-000219 

GP-303564 

Thus, the corruption of the measurement signal at higher air mass flow rates can 
be significant. 

[0026] The present invention provides a Kalman Filter (KF) -based 
signal processor. A standard KF combines the concept of recursive least 
squares estimation with a state space model and noise model to optimally 
estimate the signal value. The KF estimates the signal value based on a linear 
model of the sensor signal. The linear model of the sensor signal includes a 
linear dynamic model, a description of the covariance of the errors in the signal 
dynamics, and a description of the covariance of the errors in sensor 
measurements. The linear model is expressed in state space form as: 

Xk+1 = AkXk + BkUk (Process Model); and 
yk = CkXk (Measurement Model) 

where: 

Xk is the state vector (i.e. what is to be estimated); 
yk is the measurement vector (known); and 
Ak, Bk, and Ck are the system matrices (known). 

[0027] The KF estimates the sensor signal by using a form of feedback 
control. In short, the KF estimates the sensor signal at some time and then 
obtains feedback in the form of noisy measurements. The KF equations fall into 
two groups: time update equations (predictor) and measurement equations 
(corrector). The time update equations project the current sensor signal and 
error covariance estimates to obtain estimates for the next time step. The 
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measurement update equations incorporate new measurements into the 
estimates to obtain an improved or updated estimate. 

[0028] The KF-based signal processor of the present invention 
introduces a mechanism that propagates a mass air flow estimate (x) based on a 
3"* order model. The estimator propagates the estimate error covariance based 
on the model and how noise affects signal generation. The controller 18 
processes the KF-based signal processor and controls other components of the 
system 10 in response to the smoothed signal. 

[0029] An output-error model is used to develop an estimated airflow 
model and an estimated compressor command model. The output-error model is 
able to handle heavily corrupted measurements while maintaining a standard 
transfer function form. Although the output-error model is implemented in the 
described embodiment, it is appreciated that other linear modeling techniques 
can be used. The output-error model is an approximation of the sensor 
dynamics. 

[0030] Referring now to Figures 3A-3C, open-loop response data is 
collected from the compressor mass flow sensor system. More particularly, a 
compressor command signal is generated and mass flow rate data points are 
collected from the mass flow sensor. The compressor command signal includes 
a single step increase. The exemplary compressor command signal of Figure 3A 
includes an increase from approximately 1000 RPM to approximately 8000 RPM. 
The measured air mass flow rate from mass flow sensor 16 is illustrated in Figure 
3B. A model is generated based on a 3"^ order curve-fit of the air mass flow rate 
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data points and illustrated in Figure 3C. The air mass flow rate signal is modeled 
by the following equation: 

AirFlow = x = Q^z^ -Q2Z + Qj 
The compressor command signal is modeled by the following equation: 

CompressorCommand = u = -R^Z^ + R2Z-R3 

The coefficients of the air flow signal model and the compressor command 

signal, Q and R, respectively, may slightly vary based on the particularities of the 

compressor 14 and the mass flow rate sensor 16. The coefficients of the 

exemplary embodiment include: 

(2, =2.40197810e-6 
^2 =5.97555508^-6 
^3 =3.73093133e-6 
/?, =2.93515064 
/?2 = 2.87339029 
/?3 =093817087 

[0031] The matrices that govern the KF-based signal processing of the 

present invention are provided in state-space form as follows: 

x= Ax + Bu 
x = Cx 



A = 



/?! /?2 /?3 



1 0 0 
0 1 0 

B = [l 0 of 
C = [Q. ^3] 

The KF-based signal processing takes a previous (t = k-1) smoothed estimate 
( jc*_, ) and calculates a predicted current estimate (Jc;) at the current time (t = k). 
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A new measurement sample (3^) is used to create a smoothed current estimate 
(jc^^). The smoothed current estimate (Jc^^) is used to predict a future estimate 
{^m) (t = k+1). The amount of smoothing is determined by the Kalman gain (K), 

which is directly calculated by the measurement noise variance and the estimate 
covariance (P). In general, if the measurement device (e.g., mass flow sensor 
16) is known to be noisy, the smoothing effect will match the model result {x~) 
more closely. However, if the measurement device is not known to be noisy, the 
estimate will match the measurement (y) more closely. 

[0032] The Kalman gain is defined according to the following equation: 

where: P^^ is the noise covariance matrix and is R = 1 for this system; and 

H is the measurement matrix 
The smoothed current estimate (jc^) is calculated according to the following 
equation: 

The updated estimate covariance is defined according to the following equation: 

p;:,=ii-K,H)p- 

The update of the signal based on the model (Jc;^,) is calculated according to the 
state-space equations provided above, as: 

The covariance estimate is defined according to the following equation: 
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The initial conditions are provided as: 

P- = BQB' 

[0033] It should further be noted that the A-matrix and B-matrix are 
defined as previously discussed and the Q-matrix is the process noise 
covariance, which is set equal to 10® for the present embodiment. The Q-matrix 
represents dynamics not included in the model, which result in how well the 
estimate matches the data. Finally, ( ^ ) denotes a matrix transpose and ( a 
matrix inverse. 

[0034] Referring now to Figure 4, the steps of the KF-based signal 
filtering are illustrated. Initially, the matrices are loaded and the initial conditions 
are set in step 400. In step 402, the current estimate (Jc^) is predicted based on 

the previous smoothed estimate (x^^i). The current Kalman gain (Kk) is 
calculated in step 404. In step 406, the smoothed current estimate (Jc^^) is 
calculated based on the predicted current estimate (Jc^), the current 
measurement {y) and the Kalman gain (Kk). The future estimate (Jc;^,) (i.e., 
estimate for the next time iteration) is calculated based on the smoothed current 
estimate (jc^) and the current command (Uk) in step 408. In step 410, the 
controller 18 controls the fuel cell system 10 based on the smoothed current 
estimate (Jc^^). 
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[0035] Referring now to Figures 5A-7C, a series of graphs illustrate 
and compare the accuracy of the KF-based signal processing of the present 
invention. More particularly, Figures 5A and 5B illustrate the performance of a 
recursive Kalman filter derived according to the measurements of Figure 2 with 
data points illustrated in Figure 5A compared to the filter estimate illustrated in 
Figure 5B. Such a comparison shows off-line simulation results of the recursive 
Kalman filter. The up-transient that begins at time T=42 seconds is the step 
response used to perform system identification. It should be appreciated from 
these figures, the filter does an excellent job filtering the heavily corrupted 
measurement signal even at high flow rates. 

[0036] Figures 6A and 6B show simulated, comparative performance 
between the recursive Kalman filter derived according to the measurements of 
Figure 2 and a low-pass first order lag filter derived according to the 
measurements of Figure 2. The raw data and a plot of the filtered data for each 
filter type is show in Figure A. The filter data alone for each filter type is shown in 
Figure 6B without the data points. Input data is modeled with a low-pass 
first-order lag filter (shown as a thin line in Figure 6B) having a time constant 
equal to 0.25 seconds. Input data is modeled with the Kalman filter (shown as a 
thick line in Figure 6B). Comparison of these two lines with the data shows that 
standard low pass filter may handle noisy measurements but result in significant 
lag and impercision respective to the Kalman filter in capturing the transient 
behavior of the system. 
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[0037] Figures 7A through 7C provide a comparison for the 
effectiveness of the performance across a broad range of operating conditions of 
the recursive Kalman filter derived according to the measurements of Figure 2. 
Figure 7A shows a step change in outputs from compressor 14 with Figure 7B 
showing data from mass flow sensor 16 and the Kalman filter data overlayed. 
Figure 7C shows only the Kalman filter data. A comparison of Figures 7B and 
7C clearly illustrate that the recursive Kalman filter tracks the system output in 
the presence of heavy noise corruption of the signal from the mass flow sensor 
16 to the controller 18. 

[0038] Figure 8 provides a software diagram 800 (Simulink^'^ 
implimentation) for the recursive Kalman filter (output at Black 816) derived 
according to the measurements of Figure 2. This coding is compiled to real time 
controller code resident in controller 18 for enabling execution of the algorithm 
respective to the filter for measurements (Block 802) from the mass flow sensor 
16. Common reference between diagram 800 and the equation set forth in 
paragraph [0028] above is appreciated from the R Matrix (Block 804) of the 
measurement noise analysis and from the A, (Block 808), B, (Block 810), and H, 
(Block 812), matrices in the system identification portion of the algorithm. The 
remainder of the signal lines and mathematical operations represent a 
conversion of these equations into Simulink™ context with PT1 (Block 814) being 
the comparative low pass filter. 

[0039] The described recursive Kalman filter accurately processes 
noisy measurement signals in a manner providing clear improvement over a 
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conventional low pass filter. A number of benefits are derived from the above- 
described approach. Compact fuel cell system designs are conceptually enabled 
where mass flow sensors do not need several straight pipe diameters before the 
meter to "straighten the flow" reducing the effective turbulence. Similarly, the 
need for "honeycomb" flow straighteners at the inlet to the mass flow sensor has 
potentially been eliminated. Thus, low cost mass flow sensors are more viable 
than fuel cell systems insofar as the predictive estimation software (i.e., the 
recursive Kalman filter) accomodates turbulent flow without sacrifice transient 
response. 

[0040] Compressor pulsation is also minimized as a smooth flow signal 
is provided to control the compressor. In the fuel cell stack, durability is 
beneficially augmented insofar as air flow pulsation is minimized along with 
comensurate reduction in cackled pressure pulsation. Anode pressure thereby 
more smoothly tracks cathode pressure to consistently provide a small pressure 
drop between the cathode (air) and the anode (hydrogen) gases across the PEM. 
Such consistent pressure drop management is vital to fuel cell stack performance 
and durability. Finally, transient performance is well managed with the rapid 
broad range response enabled by the predictive estimating filter. 

[0041] The description of the invention is merely exemplary in nature 
and, thus, variations that do not depart from the gist of the invention are intended 
to be within the scope of the invention. Such variations are not to be regarded as 
a departure from the spirit and scope of the invention. 
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